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Abstract 

Friedrich's proofs for the global existence results of de Sitter-like space-times and of semi- 
global existence of Minkowski-like space-times [Comm. Math. Phys. 107, 587 (1986)] are 
re-examined and discussed, making use of the extended conformal field equations and a gauge 
based on conformal geodesies. In this gauge the location of the conformal boundary of the 
space-times is known a priori once the initial data has been prescribed. Thus it provides an 
analysis which is conceptually and calculationally simpler. 
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1 Introduction 

rn _ 

. In [HI IS] the existence and stability of vacuum de Sitter-like space-times has been discussed. 

Moreover, reference [11] provides a semi-global existence and stability result for the development 
0^ , of hyperboloidal initial data which is close to Minkowski data. These results were subsequently 

generalised to the case where the gravitational field is coupled to Maxwell and Yang-Mills fields 
in [12]. These results make use of the Einstein conformal field equations — see e.g. [8l l9l 171 HO]— 
to reformulate Cauchy problems which are global or semi-global in time into problems which are 
local in time. Given one of these local Cauchy problems, then using powerful results of the theory 
of quasi-linear symmetric hyperbolic systems, e.g. [2TJ1[H1[22] 5 it is possible to prove the existence 
of solutions which are close to some explicitly known reference solutions — the de Sitter space-time 
and the Minkowski space-time. This particular strategy to prove global and semi-global existence 
and stability only works in 4 dimensions, although it should be mentioned that alternative proofs 
[U [2] have been obtained using the so-called Fefferman-Graham conformal invariants which are 
valid in arbitrary even dimensions. 

When discussing the conformal structure of space-times it can prove valuable to make use 
of gauges based on conformal invariants. One of these invariants, the conformal geodesies, has 
been introduced in [T5] as a tool for the local analysis of the structure of conformally rescaled 
space-times. These conformal geodesies are associated to conformal structures in a similar way 
as geodesies are related to a metric. More importantly, conformal geodesies retain their character 
upon conformal rescalings. As in the case of the usual Gaussian coordinates, coordinates on a 
fiduciary spacelike hypersurface are kept constant along a fixed congruence of timclikc conformal 
geodesies. In [16] it has been shown that on the Schwarzschild-Kruskal space-time it is possible 
to construct a system of globally defined conformal Gaussian coordinates. 

In [13j a more general set of conformal field equations has been derived: the extended confor- 
mal Einstein field equations. These conformal equations are expressed using a so-called Weyl or 
conformal connection. A Weyl connection is a torsion free connection (not necessarily Levi-Civita) 
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which preserves the conformal metric and for which parallel transport preserves conformally or- 
thonormal frames and the causal nature of their vectors — that is, time- like vectors are transported 
into time-like vectors, etc. The extra freedom introduced by the use of Weyl connections allows 
one to consider gauges based on conformal geodesies in the discussion of various local and global 
issues in General Relativity. For example, in [13] this type of gauge and the evolution equations 
implied by them have been used to discuss the (local) existence of anti de Sitter-like space-times; 
in |14j conformal geodesies have been used to construct a representation of spatial infinity which 
makes possible, in principle, the detailed discussion of the structure of the gravitational field in 
this region of space-time. 

One of the properties of conformal Gaussian coordinates is that in the case of vacuum space- 
times, they provide a canonical conformal factor, written entirely in terms of quantities defined 
on a fiduciary hypersurface and the conformal time parameter. In turn, this conformal factor 
provides a priori knowledge of the structure of the conformal boundary of the space-time. In this 
article we make use of this aspect to re-examine and discuss the global and semi-global existence 
results of [11) . The possibility of following this approach is somehow implicit in the literature, 
but to our knowledge has never been made explicit. The approach discussed here is in a sense 
more natural as it exploits the full freedom contained in the extended conformal field equations 
by making use of a gauge based on conformal invariants. In the original approach used in |llj 
the conformal factor is itself an unknown satisfying a propagation equation. Consequently the 
properties of the conformal boundary of the resulting space-time have to be discussed a posteriori 
in an abstract way. For example, in the case of the development of hyperboloidal data, the 
existence of a point corresponding to future time-like infinity where the generators of null infinity 
meet is inferred from a qualitative argument involving the Raychauduri equations. In contrast, 
in the approach followed here the existence and properties of future time-like infinity follow from 
explicit calculations and the requirement that the data is close to Minkowski data. Moreover, this 
information is available without reference to the existence problem of the propagation equations. 

1.1 Structure and overview of the article 

This article is structured as follows. Section [2] discusses congruences of conformal geodesies in 
vacuum space-times and the construction of conformal Gaussian coordinates on space-times whose 
time slices are diffeomorphic to the 3-sphere, or subsets thereof. The canonical conformal factors 
associated to the congruence of conformal geodesies are also studied. Section [3] is concerned with 
the conformal de Sitter and Minkowski space-times as subsets of the Einstein cylinder and with 
congruences of conformal geodesies on these manifolds. These space-times are analysed in detail 
as they will be used as reference solutions later on. Section [4] discusses the conformal boundary 
of the development of initial data sets which are sufficiently close to Cauchy data for the de Sitter 
space-time or to hyperboloidal data for the Minkowski space-time. Section[S]is concerned with the 
propagation equations implied by the extended conformal equations and the conformal Gaussian 
coordinates. The propagation equations are expressed in terms of a space spinor formalism. In 
section [6] the de Sitter and Minkowski space-times are recast as solutions of the propagation 
equations of sectional Section [7] discusses the global and semi- global existence results implied by 
the propagation equations of section [5] These existence results are a consequence of a modified 
version of an existence and stability result for quasilinear hyperbolic systems by T. Kato. This 
modified Kato theorem was first discussed in [11) and is included here for completeness (theorem 
[1]). The existence results considered here are: global existence for de Sitter-like space-times 
when data are given either on a standard Cauchy hypersurface (theorem [2]) or on the conformal 
boundary (theorem [3]); and semi-global existence for hyperboloidal Minkowski- like data (theorem 
[3J . Section [5] provides some concluding remarks about possible generalisations of the results 
provided in the present work. An alternative discussion of the behaviour of conjugate points in 
the congruence of conformal geodesies is provided in the appendix. 
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2 Setup and gauge considerations 



Let {M-Tg^u) be a space-time satisfying the Einstein field equations with cosmological constant 

We will only be concerned with the case A < 0. Space-times such that A < and suitably close to 
the de Sitter solution will be called de Sitter-like whereas if A = and the space-time is suitably 
close to Minkowski space-time it will be called Minkowski-like. The metric g^ v is assumed to have 
signature (+, — , — , — ). 



2.1 Coordinatising manifolds diffeomorphic to § 3 

We shall work with space-times which are of the form I x S where / is an interval on R and S 
is diffeomorphic to § 3 or to a submanifold thereof. The manifold § 3 will always be thought of as 
the following submanifold of R 4 : 



(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 + (x 4 ) 4 = 1 



The restrictions of the functions x A : A — 1, 2, 3, 4 on R 4 to S 3 will again be denoted by x A . The 
vector fields 

ci = x 1 di — x A d\ + x 2 dz — x 3 &2, (la) 

C2 = x 1 d 3 -x 3 d 1 +x i d 2 -x 2 d 4 , (lb) 

c 3 = x 1 d 2 — x 2 d\ + x 3 di — x 83. (lc) 

on R 4 are tangent to § 3 . In the sequel they will always be considered as vectors on S 3 . Denote 
by dui 2 the line element obtained as pull-back of 



A=l 



to § 3 . The vector fields ci, C2, C3 constitute a globally defined frame on S 3 which is orthonormal 
with respect to dcu 2 . 

Let tf> : S — * S 3 denote the diffeomorphism connecting S and § 3 . The diffeomorphism <j> will 
be employed to pull-back the functions x A , A = 1, 2, 3, 4 on S 3 to S . Their pull-back to S 
will again be denoted by the same symbol. This system of equations has rank 4 on S, so that 
one can use a suitable choice of three of the functions x A , to obtain a coordinate system in a 
neighbourhood of a point in S. Moreover, given x A on S, one can define the vectors ci, C2, C3 as 
given by (|lap - (flc]) . In the following we will extend this frame by Co to a frame in M.. The index 
letters s — 0,1,2,3, respectively f = 1,2,3, will be specifically reserved to denote components 
with respect to this frame. 



2.2 Weyl connections and conformal geodesies 

Let 

be a conformally related metric where O is some conformal factor. Let be a smooth 1-form. 
Denote by V and V the Levi-Civita connections of g and g, respectively, and by V the Weyl 
connection for g satisfying V = V + S(b) , where 
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Then V = V + S(T) with T = G^VO and V = V + S(f) with / = b - T. Further, define 
= Qb^ = 0/ M + V M 9. The Schoutcn tensor associated to the Weyl connection V is given by 

We note that whenever the connection preserves a metric in the conformal class, i.e. it is a 
Levi-Civita connection, then L^ v j = = R[^]. 

Let efe, k = 0, . . . , 3, be a g-orthonormal frame field, i.e. satisfying g(ei,ej) = rjij, with 
rjij = diag(l, — f , — f , — f), i, j = 0, . . . , 3. Denote by and V/c the covariant derivative in the 
direction of ej~ with respect to V and V. Define the connection coefficients T/ k of V in this frame 

by V^efe = T? k ej. A conformal geodesic x(t) is obtained, together with a f-form 6(r) along the 
curve, as a solution to the system of equations 

+ S{b)^ p x x x p = 0, 

i"V„6^ - \b v S{b) x \x x = L Xft x x , 

where x denotes the tangent vector to the curve x(r) and — Ag M „ with A = 6A. In what 
follows we shall often write for ± M . Given initial data for these equations in the form x* e Ai, 
x* € T Xt M., &* e T*^M, there exists a unique conformal geodesic (x(t), 6(t)) near x* satisfying 
for given r € M 

x(t )=x*, x(t ) = x*, 6(t ) = 6*. 

Conformal geodesies are conformally invariant in the sense that if x(t) and b(r) solve the confor- 
mal geodesies equations and we define a new Weyl connection V = V+S(b), then (x(t), b(r)— b(r)) 
solve the conformal geodesic equations with V replaced by V and L^ v by L M „. In particular if 
b = b then V and V coincide and the conformal geodesic equations take the form 

x v \J v x^ = 0, 
= 0. 

2.3 Conformal Gaussian coordinates 

Let S be a spacelike hypersurface in the vacuum space-time (M,g). Let h a /3 denote the intrinsic 
3- metric of S induced by g^ v . On S choose: 

(i) a positive conformal factor 6»; 

(ii) a frame field e*fc, k = 0, . . . , 3 such that g(e*i, e*fc) = <d~ 2 rnk, 

(iii) and a I-form &*. 

Given the above initial information, there exists through each point i, £<$a unique conformal 
geodesic (x(t), &(t)) with T = r,on5 which satisfies the initial conditions x(t*) = e », &(r*) = 6*. 
These curves define a smooth-caustic free congruence in a neighbourhood W of S if all data are 
smooth. In addition, b defines a f-form on U from which one can construct a Weyl connection 
V = V + S{b). A smooth frame field and a conformal factor 9 are then obtained on U by 
solving the propagation equation 

x^V^efc = 0. 
It can be seen that g(ei,ej) = 0~ 2 i]ij on U with 

x»V l ,Q = Q{b,x), 0\ S = G», (2) 

where (•,•) denotes the contraction of a 1-form with a vector. Accordingly, the frame is or- 
thonormal for the metric g^ v = Q 2 g^ v - Coordinates x A on S can be dragged along the congruence 
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of conformal geodesies, so that if one sets x° = r, one obtains a coordinatisation of U C M.. In 
this gauge one has that 

x = v = e = d Tl T J k =0, L Qk = 0. (3) 

Such a choice of coordinates, frame field and conformal gauge will be referred to as a conformal 
Gaussian system. Let g^ = + v^v u where the pull-back of h?" is the contravariant negative 
definite intrinsic metric of the surfaces orthogonal to the congruence. The Levi-Civita connection 
of h^v will be denoted by D. Since v = eo is orthogonal to S, h^ u /j,, v = 0, . . . 3, coincides with 
the 3-metric h a p = 9»/i Q /3 a, 13 = 1,2,3. In general it will however not agree with the 3-metric 
of the surface of r = constant. 

If g is a solution to the Einstein vacuum field equations with cosmological constant A then 
for a conformal Gaussian system the conformal factor and the 1-form d k can be determined 
explicitly from equation ([5]) and the initial data. More precisely, one has that 

e*^o e = e» (\ + r(&„±») + y ^Ae; 2 + i 9 »(^,M)) , (4a) 

9* = 9 = (d»,u»)r+ \&*t 2 , (4b) 

and 

d = 6, d a = (&*, 9*e a *), a = 1,2, 3, (5) 
where we have set &* = T», r» = and used the identity 

&e = x + ^(d,d), (6) 

In the above expressions, quantities with the subscript * are regarded as constant along the 
conformal geodesies, while (/"(•, •) denotes the contravariant form of the metric g^ v applied to a 
pair of 1-forms. As long as the congruence of conformal geodesies does not degenerate, then at 
the points where 9 = 0, V^.9 7^ one finds that 

g i {d,d)=T 1 kl V k QV l Q = -2\. 

The nature of the conformal boundary (spacelike, timclike, null), defined by the conditions O = 0, 
Vfc0 ^ 0, can be deduced according to whether A is negative, positive or zero, respectively. 

2.4 Conjugate points 

As mentioned in the previous paragraphs, it is necessary to check that the conformal geodesic 
congruence does not develop any conjugate points or caustics. This would lead to a break down 
of the conformal Gaussian coordinates. For our analysis we use the conformal Jacobi fields 
rf = r] k e^, as their component hj k r\ k vanishes at a conjugate point [16,. They satisfy 



drVk = Xjkrf > ( 7a ) 
3jk + Lji 



dl m = -{E jk +L jk )rf, (7b) 



d 6 ^ = -d T {E jk )rf) + Y 0jk rf, (7c) 

where Xjk denotes the components, with respect to the frame e k , of the second fundamental 
form of the surfaces of constant r, and we have used the 4-dimcnsional Cotton- York tensor 
Yijk = V[jl/jifc, and the electric part of the Weyl tensor, Eij = Coioj- 

3 The asymptotically simple model solutions 

In this section we discuss the exact solutions of de Sitter (A = —1/2) and Minkowski (A = 0) in 
terms of conformal geodesies and initial data for the regular conformal field equations. It is well 
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known pjj] that both solutions can be conformally embedded into the Einstein cylinder R x S 3 
with metric 

g E = d 2 t-d 2 x-sin 2 xda 2 , 

where \ € [0, n] and da denotes the standard metric of § 2 . The curves x^{t) — x^+tS^ — x^+tu^ 
are conformal geodesies written in terms of the coordinate t and u — d t . Along these curves we 
have 

®e{t) = (i + ^- 
v{t) = Q e {t)- x u, 

T dQ E 
b E{T) = 7:dt= ——, 

with 

t 

t = 2 tan - , 
2 

and where @e is the conformal factor generated along the geodesic by ([2]) with Qe* = 1 at t = 0. 
Observe that r — > ±oo, &e —> oo and v — + as t — > ±7r. The metric g = & e 9e is given by 

g = Q%g B = d 2 r - (l + x) ^ x + sin ' xAa2 ) ^ 

and satisfies g(v, v) — 1. Note that all curves are orthogonal to the surfaces of constant r. In the 
rescaled space-time b — bE — O^MOs = 0, and thus, the curves are geodesies with respect to g. 

3.1 The de Sitter space-time 

The de Sitter space-time is embedded into the Einstein cylinder using the conformal factor 

'i-t^ _ 2-e E 



O n = COS t 



1 + $) ©I 



The conformal factor £Id vanishes at r = ±2, where the conformal boundary = U is 
located. Using 

G?fi£> r 

= -— = — dt, 

"15 1 — 4r J 
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setting 0£>„ = 1 and using equation we get that along the congruence of conformal geodesies 

e D ( T ) = l- Jr 2 = n D e E , 

. , ^ T -|t d(» D 9 g ) 
6 C = b E + T D = — dr = — — . 

1 - ''doe 

On the Cauchy surface r = we have &d(0) = 0. From equation (jla|) we recover Op, as above. 
We see that rescaling by Od gives the metric g once more. The Cauchy surface r = — 2 represents 
the past conformal boundary t?~ . If we redefine r — * r + 2 we get on the conformal de Sitter 
space-time that 

d D = ® D b D = --(r-2)dr, 
u 

v 



l + Kr-2) 2 ' 
with the initial data at .y~ given by 

do* — dr, (d£),u)* = l, On* = --. 



We observe that in this case On is given by the formula j4b 
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3.2 Minkowski space-time 

The Minkowski space-time will be embedded using the conformal factor 

Q M = cos(t + - ) + cosx. 

For convenience we have shifted the standard embedding by ir/2 to the past here, so that the 
usual Minkowskian hyperboloid, which is usually embedded at t = 7r/2, is now located at t = 0. 
We have 

dSl_M cos tdt + sin x dx 

Tm = -7: — = ■ 

s Im sin t — cos x 

It follows that the conformal geodesies satisfy 

t\ , cosidi + sinxdx 
5m = ob -t- T a/ = tan | - ] dt 



2 / sin t — cos % 



.,2 



COS - It, 

\2 ' 



Hence at t = = r, we get the following information on the canonical Minkowski hyperboloid 

CIm* = cosx, 
d,M* = - (dt + sin x) , 
(d,M, «)* = -1, 
/t s (d,d)* = sin 2 x- 

Substituting these into formula l]4a[> we get 



9 M = cosx 1 1- —r I = ^dQe, 

V cos x 4 / 

which vanishes at 

1 ± sin x 
t = 2 -. 

cosx 

3.3 Conjugate points in the reference solutions 

For the reference solutions discussed in this section the electric part of the Weyl tensor Ejk and 
the 4-dimensional Cotton- York tensor Yijk vanish. Furthermore, the components of the second 
fundamental form Xjfc(O) vanishes. Using equations (|7a[) - ([7cl ) one finds for the chosen congruence 

Now, 8b ^ for r £ (—00,00). Thus, the Jacobi fields will never be tangent to the curve nor 
vanish. Hence the congruence is free of conjugate points. 



4 The structure of the conformal boundary for nearby space- 
times 

In this section we use the formulae (Ha|) and (|4b|) to study the conformal boundary of space-times 
which are constructed as the development of initial data which is close to either de Sitter Cauchy 
data or to hyperboloidal Minkowski initial data. 



4.1 Space-times close to de Sitter 

As discussed in [6], for de Sitter- like space-times one can formulate two slightly different Cauchy 
initial problems: one where data is prescribed on a standard Cauchy hypersurface, and a second 
one where the data is prescribed precisely on one portion of the conformal boundary, ,y~ . 
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4.1.1 The case of de Sitter-like data away from J* 

Here we assume that we are given a space-like hypersurface S which does not intersect J^~. So 
is given by (|4aj) as 6» ^ on 5. In fact, without loss of generality we could set O* = 1. The 
conformal factor vanishes at 

-26*(d,u)* ±20*J\2\ + gt(d, d)„| 

t± = = (9) 

2X + gHd,d), 

which gives the location of the conformal boundary, J^ ± as smooth space-like hypersurfaces. If 
S is topologically § 3 — as it is being assumed here — then 

y ± = {t±} x S 

are topologically § 3 . On J !± we have Vfe9V fe 9 = — 2A and hence both hypersurfaces are space- 
like. 



4.1.2 The case of de Sitter-like data on J? 

We discuss now the data for a conformal geodesic congruence that starts on the smooth hyper- 
surface S which represents J?~ . Hence, the initial data satisfies the condition O* = on S and 
takes the form (|4bj) . Equation (J6j> implies that g\d, d)* = — 2A and thus must be time-like as 
A < 0. On the other hand, 6* is free data on S. Having set 6* = Y*, that is, = (VO)* and 
u* = n, where n is the unit normal of S with respect to g, it follows from 9„ = on 5 that 



d* = (d, v)«v* with O* = (rf, w)„ = ±y — 2A. 

We choose the positive root so that O is positive in the future of S. Thus with respect to the 
Weyl propagated orthonormal frame et we obtain from [5] 

d fe (r) = (J^2i+e,T,0,0,0j ■ 
The conformal factor vanishes at 




Hence, the location of is determined by the free data O*. On we have = — d* and 

VfcOV fe O = — 2A and hence, again, it is a space-like hypersurface. 



4.2 Hyperboloidal Minkowski data 

We discuss now how to use formula (|4ap to gain a priori information on the conformal boundary 
of the domain of dependence of hyperboloidal initial data which is close to Minkowski data. Given 
a 3-dimensional manifold S with the topology of § 3 , we consider S C S, with dS diffeomorphic 
to § 2 . Furthermore, consider a function f2 on S such that SI > in the interior of 5, and 11 = 
on 2 e dS. The function SI can be obtained naturally as part of a solution to the conformal 
Hamiltonian and momentum constraint — for a discussion and more details on this see [TUJ [TT] . 

It is noted that 

where /i B (&*,&*) < 0, and (&*,&*)] = — /i" (&*,&*) since is taken to be negative definite. In 
order to make use of formula (j4a| , the following particular choices of initial data will be made 

O* = SI, = ©*, d* a = S16* a = D a fl, 



where &* a = e^6* M , D a — e^V M = e"D a , a = 1,2,3. In particular, the sign of 0» contains the 
information about which part of the locus of points such that = should be considered as the 
conformal boundary — see the discussion below. On S\Z the function O* is determined by © 
for the initial data 0* is taken to extend smoothly to Z. On Z, the formula ([6]) implies that 
g^(d, d)» = <7"(V0, VO)* = 0. For hyperboloidal data (V0)„ ^ on Z, and thus d* must be null 
on Z. We must thus have 

S 2 = (d^v*) 2 on Z, 
where S = yj\h$(d, d)*\ = yJ\D^VtD k Vl\ on S. Consequently, 

u 4 2n 

/i p (6*,6*) = with cj = 

The functions f2, 2 and w will be extended off S by requiring that they remain constant along a 
given curve of the congruence of conformal geodesies, and will be denoted again by f2, 2, ui. 

In order to further discuss the structure of the conformal boundary, we analyse the zeros of 
0. For curves passing through Z the zeros are located at 

-0*±0* 

T± = 



0, 



whereas on S \ Z we can write 



1 n 1 

4* 



= 0,(l+«T+( 7a 2 -- 2 -)r 2 ), (10) 



where a = (b*,v*). Note that a can be chosen independently of f2 — this fact will play a role 
when discussing the existence of solutions. The function a will be extended off S in the same 
way as it was done for fi, S, oj. The roots of are given by 



T± = 



a 2 il 2 - S 2 ' 



Accordingly, one defines 

J f± = S^(T,X A ) ERXS T = T±(x A )y (11) 

This shows that the location of J^^ is predetermined by the initial data $7, d*o,-D/c^ and well 
defined as long as the congruence does not degenerate. One sees that as 0* — > one has 
t± — > 0, — 20*/0*. Thus, with the continuation of 0» onto Z described above one finds that 
and are smooth hypersurfaces, whenever dO ^ 0, and Z is the intersection of with {0} x S 
as expected for hyperboloidal initial data. In analogy to the model case of the hyperboloids in 
Minkowski space-time, the development of general hyperboloidal data has a conformal boundary 
which either corresponds to ^ + or to but not both. This information is contained in the 
sign of the free datum O*. The conformal factor is positive on the physical space-time {M.g^). 
So if 0* > on Z, then M. lies to the future of ^ . In this case one speaks of a hyperboloid 
which intersects past null infinity, and thus the conformal boundary is identified as/". Whereas 
if 0* < on Z, then M. lies to the past. Then the hyperboloid is regarded as intersecting future 
null infinity and then gives the conformal boundary. Without loss of generality, in the sequel, 
we shall only consider hyperboloids intersecting future null infinity, so that O* < on Z. Then 
^ + is given by r + as identified above. We remark that the solution t_ is of no interest to us as 
it lies outside the domain of dependence of S . 

As S is a compact set, there is a point in the interior of iS for which D^l = and hence 
5 = 0. If the data is close enough to Minkowski data this critical point is unique. If one makes 
the choice a = 0, then one finds that 

t± = ±— ► oo as S — > while £1 > 0. 
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Thus, in order to have a conformal representation of the domain of dependence of S for which r 
remains overall finite on J? one needs a^O. 

In order to identify points which can be regarded as representing timelike infinity, one needs 
to investigate the critical points of O, that is, the points for which dO = 0. One has that 

<1B = ( I +,.r+ ( - { o 2 - ^ J - J ) dfl 

-ft I a + 2t I i« 2 ^ ) ) dr + ttrda + Sir 2 ( -ada + -^-dw 



In particular, we are interested in analysing the conditions d0 = on J^ ± , with fi ^ 0. By 
construction one finds 







to 




^ 2 , 





0. 

A necessary condition for the vanishing of d0 on J^ r± for ^ is that 

a+2T± (r 2 -i)=°- 

A short calculation shows that the latter is equivalent to 

S 2 = D k VLD k VL = 0, 

so that D k Vl = 0. Now, if D k fl = then r± = — 2/a. Note that r± > if a < — that is, if 
O* < 0. Thus, in order to consider a conformal representation which includes the point i + , one 
needs to consider a^0. This condition will be assumed in the sequel. Let 

T~i+ = -2/a. 

From the discussion in section 13.21 it follows that in particular for the development of Minkowski 
data one has that Tj+ = 2. 

Another computation shows that 

VLr i+ da + nrf+ Qada + ~^ dw ) = 

if dfi = 0. To show this, one uses the fact that if D k £l — and £1 ^ then it follows that 
Dfc(S 2 ) = and furthermore, that 1/ui = S/2J7 = 0. In view of this, one defines i + G M x S as 
the unique point for which r = T t + and dQ = 0. 

To conclude the discussion of the point i + , we look at the Hessian of O at i + , using the 
conformal Gaussian coordinates (r, x A ). Using D k Q = and D k {E?) = from above, one has 



V -4 V » ( - ) = 40 2 at * 



Thus wc find that 



V^V B 6| l+ = fiV^Vg ( ar + ( -a 2 - \ ) r 2 



W 2 



1 1 V V " 2 
<b V 4 Vkm fir 2 ( -V..,<»V b -.» • ^V. 4 V (; n 



Consequently 

2 

V7 „Y7„<=> — _ 



V^V B 6 = -^V^aVfja - ^ I ■ on i+ - ( 12 ) 
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Similar direct calculations render 

V.4V06 = V V^6 = -nV A a, on i + , (13a) 
V V 6=i^a 2 , oni+. (13b) 

Thus, if one chooses the function a such that a ^ (in order to have finite) and such that 
the 4x4 matrix with entries given by equations (p~2|) . (|13a|) and (|13b|) has no zero eigenvalues, it 
follows that the Hessian of on i + is non-degenerate. Accordingly, the point i + can be rightfully 
regarded as the timelike infinity of the development of the hyperboloidal data prescribed on S. 



5 Evolution equations 

The general conformal Einstein field equations introduced in [T5] — see also [13 E] — are a 
generalisation of the original conformal equations which allows for the use of Weyl connections. 
The use of Weyl connections makes it possible to consider more general gauges when one is 
confronted with the need of deriving a system of propagation equations out of the conformal field 
equations. In particular, it makes it possible to make use of the conformal Gaussian coordinates 
discussed in the previous sections. 



5.1 Evolution equations in a frame formalism 

As shown in, for example [17j , the general conformal field equations together with the gauge given 
by J3J imply the following system of propagation equations: 

8 e s — -f p e s 

= ~ f ^ f /o + - S\t i0 - g l0 L 3 k + Qd 0j \, 

d T Lji = —Tj P Lpi + o?fcd 0j fe j, 

v fe ^ fe ; = 0, 

where dijki = 0~ 1 G\jki denotes the components of the rescaled Weyl tensor with respect to the 
frame eg, while the conformal factor O and the 1-form dk are given by (|4a]) -(|5j). 



5.2 Evolution equations in a spinor formalism 

In view of future applications, instead of the frame evolution equations discussed in the previous 
section, we shall use a spinorial version thereof. Using a spin dyad {o, l} such that 

T AA ' =V2v AA ' = o A o A ' + l A I A ' , 

the evolution equations take the form 

fox » — r EE ' h- dd ' c s 

V ZU T e qqi — — 1 cc , DD' T e EE'i 

\F)H V A — -ft F f A -4-f F f A ~)t CC ' 

+l dd , bc >t ac ' + ecj) A BCD T c D ,, 

\/2d T L AA , cc , = —(t AA , F B L FB , cc , + T AA F b ,L bf , cc ,)t bb 

— d EE (4 , EABC^E'C'T B A/ + 4>E>A'B>C'£ECT A B ), 
V2d T <j)ABCD = T FF V F'(d4>ABC)F ~ l"F'(flV FF 4>ABC)F, 
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where one has the correspondences (via the Infeld-van der Waerden symbols) 

e ' \ l— * e ' CC'i 

f fc f AA' _ f A A' , f A' A 
S'i X BB' CC SB' C e C" ' SB' C' e C i 

Lij i ► Laa'bb 1 -, 

dijkl l-> dAA'BB'CC'DD' = ^ABCD^A'B'^C'D' + ^A'B'C'D'^AB^CD > 
rfj l-> ^AA' , 

and the factors of \/2 arise from the normalisation taa'T' 4 ' 4 = 2. 
5.3 Evolution equations in space-spinor formalism 

Next, one introduces a space-spinor formalism by using the spinor to eliminate all primed 
indices and then one splits the equations into symmetric and skew parts. 

A space spinor Qabcd = ©ab(cb-) 1S introduced and then decomposed such that 

&ABCD = L C C'AA'T A B T° D = Q(AB)(CD) + ^AB^Q (CD)- 

For the spin coefficients Taa'bc, one observes that Taa'bc = r^A'sc + ^abJca 1 and defines 

Taboo = t b b Tab'CD, 
which in turn, will be decomposed as 

TaBCB = (£,ABCD ~ X(AB)Cd) ~ ^AB.fcD- 

The spinors in the latter equation possess the following symmetries 

TaBCB = ^AB(CD), XABCD = XAB(CD), S.ABCD = £,(AB)(CD)- 

The term £abc d is related to the intrinsic connection of the leaves of the foliation defined by taa 1 ; 
the term X(AB)cd to the second fundamental form of the leaves; and Jab to the acceleration of 
the foliation. The spinors £abcd, X(ab)cd and Jab are calculated from Tabcd via the relations 

S.ABCD = ABC D + T B B T c ° T D D TaB'C'D'^ , 

XABCD = ( T B B T C' C T D D TaB'C'B' - ^ABCD^j , 

JAB = -T 1 CC'AB- 

The frame fields e s AA , are decomposed using 

e AA' = ~^ TAA ' ~ rB A' e AB: 

r _ Br 
e AA' — ~ T A' e AB> 

with 

e AB = T (A e B)B'- 

The fields e s AB are associated to spatial vectors, and hence, they satisfy the reality conditions 

e A B = -r A A 'r B B 'e% B ,. (14) 
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Using the gauge given by ((3|) it can be shown that the extended conformal field equations given 
in [14] imply the following evolution equations for the unknowns e AB , £,abcd, /ab, X(AB)cd, 

Q(AB)CD, @G G CD : 

dre Q AB = -~X(AB) EFe EF - Jab, (15a) 
dre r AB = -X(AB) EFe EF, (15b) 

d T ^ABCD = ~X(AB) Ei ? ^EFCD + ^J=(^ACX(BD)EF + ^BDX{AC)EF)f EF 



-V2X(AB)(C E fD)E - y{ e AC®F F BD + e BD®F F Ac) ~ ^VABCD, (15c) 



drlAB = -X {A B) EF fEF + -^=Q p F A B> ( 15d ) 



EF 

9tX(ab)cd = ~X(ab) Xefcd - Q(cd)ab + Qvabcd, (15e) 

d T Q(AB)CD = ~X(cd) EF ®(ab)ef — d T <driABCD + iV2d E \ a ^b)cde, (15f) 

dr^AB = X(AB) EF EF + V2d EF f ]A BEF, (15g) 



where 



VABCD = 2 (^ABCD + T A A Tg B T C C T D ° <j> A > B'C D')i 
HABCD = --^{4>ABCD ~ T A A T B B T C C T D D <j> A >B'C>D>)i 

denote, respectively, the electric and magnetic parts of of 4>abcd- 

The evolution equations for the spinor 4>abcd are derived from the Bianchi equations. De- 
pending on the need, several alternative systems can be deduced. Here, we will consider the one 
which was called the standard system in references [131 114) . Let 



\BCD — foe ABC D + 4>l e ABCD + ^I^ABCD + <l>3 e ABCD ~ 4>4<^ABCD^ 



where 



k = <t>{ABCD) v i = 0,...,4, 



_ 1 _ 2 _ 

e ABCD = °(AObOcOd), e ABCD = l (AObOcOD), e ABCD = l (A>-BOcOD), 

3 4 

e ABCD = 1 (A l BI<COd), ^ABCD = i (A^B^C^D)- 

In the previous expressions the subindex (ABCD)i indicates that after the symmetrisation i 
indices are set to zero. One has the following Bianchi propagation equations: 

(V2 - 2eQ 1 )<9 r o + 2eo <9 T </>i - 2el 1 df<t>Q + 2eo <9r</>i 

= (2r oii - 8rioio)0o + (4r 00 oi + 8T WO a)(t>i - 6r OO oo02, (16a) 
V~2d T (j)i - e? x 5 T ^o + eo(A02 - e^df^o + e r 00 d f ^2 

= — (4Txno + /u)</>o + (2r 00 ii + 4rnoo - 2foi)<t>i + 3/00^2 - 2r OO oo03, (16b) 
V2d r (f)2 - e\ r d T <l>i + e° o <9 T 3 - e^Sf^i + e^drfa 

= -Tun4>o - 2(r 110 i + fu)(f>i + 3(r on + rnoo)02 

-2(r ooi - /oo)<^3 - r O ooo04, (16c) 
\/2<9 T 3 - e? x 5 T ^2 + eo O 9 T 04 - e^df^z + e r 00 df<t>i 

= -2r 1111 1 - 3/ n 02 + (2r U00 + 4r 0011 + 2/ ol )0 3 - (4r 0001 - / oo )0 4 , (16d) 
(V2 + 2e% 1 )d T <f> 4 - 2e° n d T <f> 3 + 2e F 01 dr<p 4 ~ 2ef 1 $.fo 

= -BTuufa + (4r mo + 8r om )0 3 + (2r U00 - 8T QW1 )(f> 4 . (16e) 
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5.4 Evolution equations for the Jacobi field 

To the above propagation equations we will have to append an evolution equation for the Jacobi 
field. The conformal Jacobi field has a spinorial counterpart tjaa 1 which can be split as 



VAA> 



1 



-VTAA' ~ T A'VAB, 



with 



f] = TjAA'T 



_AA' 



1]AB 



r (A 7 lB)B>- 



Conjugate points in the congruence of conformal geodesies arise if t]ab = 0. The components rj, 
rjAB satisfy the propagation equations 



V2d T rj = f A BV 



AB 



V2d T r/AB = Xcd(AB) 7 T 
5.5 Structural properties of the evolution equations 



(17a) 
(17b) 



We discuss now some general structural properties of the equations (|15a[) -( fT5*g| ), (|16b[) - (|16e[) . 
(|17ap - (|17b[) which will be used systematically in the sequel. Introduce the notation 

v = {e ABl T ABCD-, & ABC d) , = (00,01,02,03,04) , 

where it is understood that v contains only the independent components of the respective spinor 
— which are obtained by writing linear combinations of irreducible spinors, as discussed in the 
previous section. The unknown vector v has 45 independent complex components, while (f> has 5 
complex components. In terms of v and 0, the propagation equations (|15ap - fT5"g| and (|17ap - (|17b[) 
can be written as: 

d T v = Kv + Q(v, v) + L(f>, (18) 

where K and Q denote , respectively, a linear constant matrix-valued function and a bilinear 
vector-valued function both with constant entries and L is a linear matrix-valued function with 
coefficients depending on the coordinates. Similarly the system (|16b[) - (|16ep can be written as 



V2£<9 T + A AB e AB d s 4> = B(T 



ABCD) 



(19) 



where E denotes the 5x5 identity matrix and A AB ( 



AB> 



s — 0, ... ,3, are 5x5 matrices depending 



on the coordinates, while B(T abcd) denotes a constant matrix-valued linear function of the 
connection coefficients V abcd- For later reference it is noted that 



/ V2 



V2E- 



A ab p° 

A e AB 



ze 01 


zc 00 











e ll 


V2 


pO 
e 00 











-e° n 


V2 


pO 















pO 
e 00 











-Hi 


V2 + 2e 



01 



and that 



A AB e AB 



with r = 1,2, 3. From the reality condition (fl4"]) one has that 





2 e oi 


2 e oo 











\ 






V2 


e oo 















-ell 


V2 


e 00 

















V2 


p f 

L 00 




V 











-2e\ x 


Hi 


) 



p° 

c 00 



-01 



e ll, 



-p r 
e ll, 



-° 



01' 



-01 



-01' 
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so that in particular and are real. Strictly speaking, a discussion of the symmetric 
hyperbolicity of the system (fT8|) -(fT9 f should be carried out using real unknowns. In order to ease 
the presentation in the sequel, we write 



ib s , 



-01 



where a s , b s and c s denote the components of real vectors. Thus, making use of the splitting 
4>j = Rc(0j) +ilm((j)j) and by multiplying the equations (|16a|) - (|16e|) by suitable numeric constants 
one finds that the 5x5 matrix \/2E + A AB e° AB implies the 10 x 10 matrix: 



V2 



I°(a ,6°,c ) 



V 



a 







-b° 







V2 
a 




b° 


-b° 








a 
V2 
a 





6° 


-b° 








a 
V2 

a° 






b () 


-6° 



1 
71 





6° 












b° 



1 

71 



ft" 








b° 





V2 
o° 








-6° 

6° 




a 
V2 
a° 








-6° 


6° 




a 
V2 








-b° 










1 

71 



In particular note that if one sets a = b° = c° = 0, then one gets 

i°(0, 0, 0) = diag V2, V2,V2,-L i, ^2, ^2, V2, . (20) 

Similarly, from the 5x5 matrix A AB e r AB one deduces the (real) 10 x 10 symmetric matrix: 



/ 


-e 


a f 














-V 











or 


V2 


a f 








// 





-¥ 













a f 


V2 


af 








b f 





-b f 













af 


V2 


a f 








& f 


















af 


c f 











1/ 










V 











-e 


af 













-b f 





b f 








af 


V2 


af 













-b f 





b T 








af 


V2 


a f 













-b f 





6 f 








a f 


V2 


af 


V 











-b f 














of 


c f 



A f {a\b\c f ) 



) 

For each s = 0, 1, 2, 3, the matrices A s (z) have entries which are polynomials of at most degree 
one in z = (a s , 6 s , c s ). We can rewrite them in the following form 

A s (z) = A s (0) + A s (z) 

where A s {0) = i s (0,0,0) and A s {x + y) = A s {x) +A s (y). 



6 The conformal de Sitter and Minkowski space-times as 
solutions to the conformal field equations 

In the conformal geodesic gauge given by ([3]) both the de Sitter and the Minkowski space-times 
are conformally rescaled to the unphysical space-time (M,g) where g = Q%gE — see equation 
([8]) . The connection and curvature spinor components form the variables v and <j) that satisfy the 
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equations (|18|) and (fl9|) — respectively (|15aj) -( |T5g| ) and ()16a,jl - (jl6e[) — can be directly calculated 
from the components for the Einstein cylinder. More precisely, a straightforward calculation using 
the results of section [3] renders 





(21a) 


f 4 f 

e AB - T 2 +4 Cr AB> 


(21b) 


Sab = 0, 


(21c) 




(21d) 


£,ABCD — 9 , hABCD, 
T A + 4 


X(AS)CI5 = t 2 , ^h-ABCD, 


(21e) 


-2 

6,4500= , , .hABCD, 
T -\- 4 


(21f) 


4>ABCD = 0, 


(21g) 



A' 

where tr^ B = ^mT b j are the spatial Infeld-van der Waerden symbols and 

hABCD = 

The solutions to the conformal Jacobi equations (|17ap and (|17b|) for the reference solution are 
given by 

?/ = 0, (22a) 
^B=(l + ^W B . (22b) 

It is important to note that expressions (|21a[) - pTg| and (|22a|) - (|22b)) are valid for both the 
conformal de Sitter and the conformal Minkowski space-times. What distinguishes these two 
conformal space-times is the use of the appropriate conformal factor Qd, resp. @m — cfr. section 
El It should also be observed that the expressions ()21a|l -( [2Tgl ) and (|22a| - (l22bj) are analytic 
functions of r for r G M. Important in the sequel will be that rf AB are non-vanishing for t € 
(—oo, oo). 

7 The existence and stability results 

Having written the evolution equations in the conformal geodesic gauge, we proceed now to 
analyse the existence of solutions close to the explicit reference solutions of section [6] Following 
the original discussion in [11] , the desired existence and stability results are obtained making use 
of slight modifications of very general theorems by Kato on the properties of symmetric hyperbolic 
systems [217 1 [21 1 [22] . 

7.1 Some further structural properties of the propagation system 

Let u = (Re(f), Im(?j), Re(</3), Im(</>)) with v and <f> as in section [531 The unknown u takes values 
in M N for some N S N. The evolution equations (jl5ap - fT5g| ) and (|16a|) - (|16e|) — or their matricial 
counterparts (fT51) and (fH?|) — render a system of quasilinear partial differential equations for u 
which has the form 

3 

A°(u) ■ d T u + Af i u ) ■ c f( u ) + B ( r > ^ u) • u = 0, (23) 
f=i 

with Cf(u) denoting the vector fields (|lap - (jlcp acting on the unknown u. Furthermore, 

■ 4 °<«> - ( o i»°(„) ) • • 4 '<»> - ( o !'(„) ) • 
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with E, denoting, respectively, identity and zero matrices of the suitable dimensions and A°(u) 
and A r {u) as given in section 1531 In particular given any z £ M. N , the matrix valued functions 
A s {z), s — 0,1,2,3 have entries which are polynomial in z. These polynomials are at most 
of degree one and have constant coefficients. The matrices are symmetric t (A s (z)) — A s (z), 
z € M. N . The matrix valued function B = B(t,x a , z) with (r, x A , z) 6 M x S x M N , has entries 
which are polynomials in z with coefficients which are analytic functions on K x S. Note that 
as S is diffeomorphic to § 3 , then B can be regarded as a matrix valued function with domain 
R x § 3 x R^ — this point of view will be adopted in what follows. These polynomials are at most 
of degree one. Following the decomposition in section [5.51 one can write 

A*(z) = A*{0) + A\z), 

B(t, x a , z) = B(t, x a , 0) + B(t, x a ,z), 

with 

A s (y + z) = A°(y) + A s (z), 

B(t, x a , y + z)= B(t, x a , y) + B(t, x a , z). 

Let u' denote the explicit reference solution given by (|21a[) -( f21g[ ). Set 

u = u' + w. (24) 

This is in essence a definition for the new unknown w — which gives the perturbation from the 
reference solutions. Substitution of the Ansatz (|24|) into the system (|23| yields 

3 

A°(u' + w) ■ d T w + Af ( u ' + w ) ■ c f( w ) + B ( r ) x A ,u' + w) -w + A a (w) ■ d T u' = 0, 
where it has been used that 

3 

A (it') ' d T + ' c ^ u ') + B ( T > xA ' «')■« = °> 

and that Cf(u') = — the reference solutions have no spatial dependence in the gauge being used. 
Thus, w satisfies 

3 

A°(u' + w)-d T w + ^2 A f {u' + w) ■ c f (w) + B(t, x a , w)-w = 0, (25) 
r=i 

with B(t,x a ,z) again a matrix valued function with entries which are polynomials of at most 
degree one and coefficients which are analytic functions on 1 x S 3 , such that 

B(t, x a , w)-w = B(t, x a , u +w)-w + A°(w) ■ d T u + B(t, x a , w) ■ u . 

We define new matrices A°(w) = A a (u' + w) and A r (w) = A r (u' + w) so that 

3 

A°(w) + ^(^) • cr(w) + B(t, x A , w) ■ w = 0, (26) 
f=l 

For later use it is noted that the matrix A°(u') = _4°(0) is diagonal with constant entries. More 
precisely, it has the form 

- ( i"(0) ) ■ 

with E an identity matrix of the appropriate dimensions and A (0) as given by equation (|20[) . In 
particular, all the entries of A°(u') are bigger or equal to l/y/2. 
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7.2 The Kato existence and stability result 

Let D and d/i denote, respectively, the Levi-Civita covariant derivative and the volume element 
associated with the standard metric on § 3 and D f denotes the covariant derivative in the direction 
of Cf. On the space C°°(g 3 ,R N ) of smooth Revalued functions on § 3 define for m G N the 
Sobolev-like norm 




The notation 

\D°w\ 2 = \w\ 2 , \D k w\ 2 = J2 \D fl ---D fh w\\ 

fi,...,ffc=i 

is used, with | • | the standard Euclidean norm on R N . Given to G N let H m (S 3 ,R N ) be the Hilbert 
space obtained as the completion of the space C°°(S 3 , R*) in the norm ([27]) . The unknown 
w = w(t, x) will be regarded as a function of r which takes values in H m (S 3 ,R N ). For S € M, 
to G N with < <5 < l/y/2, m > 2, set 

Df = {w£H m (S 3 ,R N ) | (z,^°(w)z) ><5(z,z),VzeR JV }, 

where (•, •) is the standard scalar product on R . Important for our purposes is that it contains 
a neighbourhood of the origin of H m (S 3 ,R N ) as the entries of ^4° (0) are bounded from below by 
1/V2. 

The original existence and stability results by Kato for symmetric hyperbolic systems of the 
form ([26]) which can be found in [22] — see also [20] [21] — have been worked out for the case 
when the frame fields c f commute. In the case under discussion where the underlying leaves of 
the foliation of the space-time have the topology of S 3 , the frame fields c f do not commute. As 
discussed in [11], Kato's result can be modified to handle frame fields whose commutators are 
those of 0(3). For completeness we quote the result given in [TT] . 

Theorem 1. Suppose to > 4, D is a bounded open subset of H m {S 3 ,R N ) with D C Df . If 
wq G D is given as initial condition for the system \2b}) , then: 

(i) There exists a T > and a unique solution w(t) of equation 126\) defined on [0, T] with 
u>(0) = wq and 

w G C{0,T;D)nC 1 (0,T;H m -\S 3 ,R N )). 

(ii) There is an e > such that one value for T can be chosen common to all initial conditions 
in the open ball B s (wq) with centre wq and radius e, and such that B e (wo) C D. 

(Hi) If the solution w(t) in (i) exists on [0, To] for some Tq > 0, then the solutions to all initial 
conditions in B e {wo) exist on [0,To] if e > is sufficiently small. 

(iv) If e and T are chosen as in (ii) and Wq G B £ (wq) with \\ Wq — wo \\ m ^ as n — > oo, then 
for the solutions w n (r) with w n (0) = Wq it holds that \\ w n (r) — w(t) \\ m — > uniformly in 
t G [0, T] as n — > oo. 

Remark 1. The point (i) in the theorem establishes the local existence of solutions to the 
equation (|26|) for sufficiently small data (but not exclusively). By (ii) there is a non- vanishing 
existence time common to all solutions arising from data in a suitably small neighbourhood of 
the 0-data — the reference solution. In particular, by (iii), if a reference solution is known to have 
a certain existence time To, then all solutions arising from data sufficiently close to the data of 
the reference solution have the same existence time. Finally point (iv) states that data close to 
a certain reference data gives rise to developments which are also close to the reference solution 
— i.e. stability. 

Remark 2. A direct computation shows that on [0, T] the solution to ([26]) is of class 
H m ((0,T) x § 3 ) C C m ~ 2 ([0,T] x § 3 ). The convergence stated in (iv) is uniform on [0, T] x § 3 . 
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7.3 Application to de Sitter-like space-times 



7.3.1 Standard Cauchy data for de Sitter-like space-times 

We start by considering the case where on a standard (spacelike) Cauchy hypersurface S one is 
given de Sitter-like initial data. The initial data consists of the values of the spinorial fields e AB , 
/ab, £,abcd, X(AB)CDi Qabcd and 4>abcd on the initial hypersurface S. This information can 
be calculated from a solution to the conformal Hamiltonian and momentum constraints. 

If the data are close to de Sitter data, then on 5 (i.e. r = 0) one has that 

e° AB = (28a) 
e AB = a AB , (28b) 
!ab = 0, (28c) 

£,ABCD = —-jhABCD + t,ABCD, (28d) 
X(AB)CD = X{AB)CD, (28e) 

Qabcd = —-jh AB cD + Qabcd, (28f ) 

4>ABCD = &ABCD, (28g) 

where quantities with a " denote quantities which vanish for exactly de Sitter data. 

Following the discussion in sections 13 . 1 1 and [4.1 . 11 initial data for the congruence of conformal 
geodesies will be chosen such that 6* = 1, b* = 0, and hence 9* = {d 7 v)* = 0. Upon this 
choice of data for the congruence, one has that the location of the conformal boundary of the 
development is given by ((9]) to be r — ±2 as in the case of the de Sitter space-time. Finally, the 
data (|28a|) -( |2"gg| ) is supplemented by data for the Jacobi field tjaa' of the form 

r? = 0, rf AB =a AB . (29) 



From the previous discussion and theorem [T] one has the following existence and stability 
result. 

Theorem 2. Suppose m > 4. Let ua = u' + uq be standard de Sitter-like Cauchy initial data. 
There exists e > such that if \\ uq ||m< £ then there is a unique solution v! + u to the conformal 
propagation equations A15a\) - [T5g\ ! and U6a]) - D^e\) with minimal existence interval r G [— 2,2] 
with u G C m ~ 2 ([— 2,2] x S) and such that the associated congruence of conformal geodesies 
contains no conjugate points in [—2,2]. The fields u' + v, imply a C m ~ 2 solution to the vacuum 
Einstein field equations with positive cosmological constant for which the sets = {±2} x S 
represent future and past conformal infinity. 

Proof. Local existence to the system of the form (|26j) implied by the propagation equations 
(|15ap -( fT"5g| ), (|16ap - (|16e[) and (|17a[) - (|17b|) follows from point (i) in theorem Q] and the observation 
that if e > is suitably small then u G D™. The reference solution given by (|21a[) -( f21g[ ) 
has existence interval (— oo,oo) D [—2,2], and the Jacobi fields associated to the congruence of 
conformal geodesies never vanish. Thus, from (ii) and (iii) in theorem [T] one has that for suitably 
small e > the developments of all data such that || u \\ m < e have a minimum existence interval 
[—2,2]. By reducing e, if necessary, one can ensure that tjab ^ for r G [—2,2] so that no 
conjugate points arise. □ 

7.3.2 Asymptotic Cauchy data for de Sitter-like space-times 

In the case of asymptotic Cauchy data for de Sitter-like space-times, whereby information is 
prescribed on an initial hypersurface S which will be regarded as the past conformal infinity of 
the development, one has that the initial value of the fields e r AB , Jab, £abcd, X(ab)CDj Qabcd 
is calculated from the following fields on a 3-metric h a p, X a a an( ^ a symmetric trace-free 
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tensor field, d a $, satisfying D a d a p = 0, where D denotes the Levi-Civita connection of h a p. The 
form of the data for the propagation equations is formally identical to that of the data (|28ap -( [28g| . 

Following the discussion of sections 13.11 and 14.1.21 initial data for the congruence of conformal 
geodesies is chosen, without loss of generality, such that (d, «)* = 1 and 0* = —1/2. The data 
for the Jacobi field is taken to be identical to that in equation (|29|) . 

The corresponding existence, uniqueness and stability result for this case is the following. 

Theorem 3. Suppose m > 4. Let uq = u' Q + uq be asymptotic de Sitter-like Cauchy initial data. 
There exists e > such that if \\ uq \\ m < £ then there is a unique solution u' + u to the conformal 
propagation equations M5a\) - (T5g\ ) and il6a]) - n~6e}) with minimal existence interval r G [0,4] with 
u G C m ~ 2 ([0,4] x S) and such that the associated congruence of conformal geodesies contains no 
conjugate points in [0,4]. The fields u' + u imply a C m ~ 2 solution to the vacuum Einstein field 
equations with positive cosmological constant for which the sets J*~ = {0} x S and .y + = {4} x 5 
represent, respectively, past and future conformal infinity. 

Proof. The proof is identical to that of theorem [2j 
7.4 Application to Minkowski-like space-times 

In the case of hyperboloidal initial data one starts with a solution (O, S, h a p, Xaff) to the A = 
conformal Hamiltonian and momentum constraints: 

2flD a D a fl - 3D a nD a n + ^n 2 r - 3£ 2 - ^il 2 ( X 2 - X a pX a0 ) + ^^X = 0, 

n 3 D a {n- 2 Xaf }) - n (d 0X - 2ir 1 £> /3 £) = o, 

where r denotes the Ricci scalar of the metric h a p and X — h a ^Xap- The solution to the conformal 
Hamiltonian and momentum constraints satisfies 

Vt = 0, S < 0, h^Dfl, DQ) = -E 2 on Z = OS, 

fl>0 on s\z. 

The value of the fields £,abcd, X{ab)cd are calculated from (O, S, h a p,Xa/3)- We set 0* = CI, 
= £ and e r AB , f ab as outlined before. The initial values of Qabcd and <f>ABCD on the initial 
hyperboloid S are then calculated from the above conformal constraint equations. 

Theorem [T] gives existence for symmetric hyperbolic systems of the form (|26|) with data pre- 
scribed on an initial manifold, S, which is topologically S 3 . Consequently the data on S has 
to be extended to data on the whole of S. Noting that S is diffeomorphic to S 3 , we consider 
in what follows 5 as a subset of S 3 . As discussed in [TT], there is a linear extension operator 
E : H m (S 1 R N ) — > H m (S 3 ,R N ) such that if v G H m (S,R N ) then (Ev)(x) = v(x) almost every- 
where in S and || Ev || m < K || v || -5, with K a constant which is universal for fixed m. As in the 
cases of the de Sitter- like space-times, the data for the equation l[26|) takes the form uq = u' Q + ua. 
The vector u' is defined as in equations ((28aj)-([28gf and thus it is defined on the whole of S. On 
the other hand, the vector uo is only defined on S, and then needs to be extended. We define the 
extended data uq by 

u = u' Q + Eu . 

It should be mentioned that the extension of the data is, in principle, non-unique and that in 
general uo will not satisfy the conformal constraint equations on S \ S. This will not have an 
effect on the development of the hyperboloidal data as D + (S) n I + (S \ S) = 0. 

In the case of Minkowski-like data we shall make use of a conformal factor O of the form given 
by formula (fTU|) with a < 0. For a sufficiently small ball of data centred on Minkowski data this 
implies that the time Tj+ is close to the one of Minkowski, namely r + = 2. The spatial location 

of i + is given by the condition D^D, = and will be close to that of i\ f . 

The existence and stability result for Minkowski-like hyperboloidal data is the following. 
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Theorem 4. Suppose m > 4. Let uo = Mq + uo 6e Minkowski-like initial data. Given Tq > 2 
£/iere exists e > suc/i i/iai 

(%) /or || Mo ||m< £ i/iera i/iere is a solution v! + u to the conformal propagation equations (15a]) - 
J_?5g[ ) cmc? (16a]) - [TBe]) with minimal existence interval r £ [0, To] and u £ C m_2 ([0, To] x S); 

(ii) the associated congruence of conformal geodesies contains no conjugate points in [0, To]; 

(Hi) for every uq £ B £ (0) there is a unique point in S\Z such that Dk^l = 0; 

(iv) for all Uq £ B £ (0) we have Tj+ £ [0,T ]. 

The solution u' + u is unique on D + (S), the domain of dependence of S and implies a C m 
solution to the vacuum Einstein field equations with vanishing cosmological constant for which the 
set ,y + as given by (11]) represents future null infinity while the point i + given by the conditions 
t = and Dk^L = represents timelike infinity. 

Proof. As before, local existence to the system of the form ([26} implied by the propagation 
equations (|15aP -( fT"5gl ) and (|16a|l - (116c|l follows from point (i) in theorem [1] and the observation 
that if e > is suitably small then u' + Euq £ D™. The reference solution given by (|21a| -( [2Tg| 
has existence interval (— oo,oo) D [0, To] D [0,—2/a], and the Jacobi fields associated to the 
congruence of conformal geodesies never vanish. Note that To is chosen independently of the 
initial data. Thus, from (ii) and (iii) in theorem Q] one has that for suitably small e > the 
developments of all extended data such that || Eu$ ||m< &/K have a minimum existence interval 
[0, To]. By reducing e, if necessary, one can ensure that ijab / for t e [0,T ] so that no 
conjugate points arise on S x [0, To]. In particular, this also holds for V + (S) C S x [0, To]. The 
characterisation of the conformal boundary then follows from the discussion of section 14.21 □ 

8 Conclusion and remarks 

In this article we reinvestigated the problem of de Sitter-like and Minkowski-like space-times with 
the help of their conformal structure. The use of conformal Gaussian coordinates has several 
advantages. Their construction is conformally invariant and provides a natural conformal factor 
along the congruence, which for vacuum can be calculated explicitly. Furthermore, the associated 
Weyl connection yields a gauge in which the extended conformal field equations are simplified 
and a symmetric hyperbolic system is obtained. It was shown that for these vacuum space-times 
the location of the conformal boundary J> + U i + can be read off directly from the initial data. It 
should be mentioned that this formulation of the initial value problem for the conformal Einstein 
field equations is amenable to numerical implementations, and indeed, a frame version thereof 
has been used in the numerical investigations of cosmological space-times described in [31 H] . 

For our calculations certain choices of initial data were motivated by the behaviour of the 
exact solution and their conformal embedding into the Einstein cylinder. Hence on the initial 
surface we set bk = It should be noted that not all congruences of conformal geodesies are 
suitable as they can develop conjugate points before or at the conformal boundary. A simple 
example is provided by the standard time-like geodesies in Minkowski space, which are conformal 
geodesies with bk = 0. In the conformally compactified picture they converge at i + , where the 
congruence has a conformal conjugate point. This is the only point of the conformal boundary 
that is reached. Moreover, it takes infinite time r to get there and the induced conformal factor 
is constant along each curve. 

Other choices of initial data are related to the parametrisation of the congruence, e.g. for 
de Sitter-like space-times 8* is free datum on . These choices affect the location of in 
our conformal Gaussian coordinates and as long as we avoid r — > oo before reaching we can 
specify them suitably. However, these choices do not reduce the class of de Sitter-like space-times 
for which theorem [3j holds. 

Finally, it is mentioned that the methods discussed here can be adapted to discuss the stability 
of other suitable reference solutions like the purely radiative space-times of The detailed 

discussion of this generalisation will be presented elsewhere. 
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A An alternative discussion of the conjugate points 

In the main part of the article the variable u = (v,<f>) contained the Jacobi fields and thus 
guaranteed that one can avoid conjugate points. In this appendix we outline an alternative 
approach. Kato's theorem — cfr. theorem [T] — can be used for an unknown u which does not 
contain r\AB- It states that one can find initial data close to the original data such that the 
components of the curvature tensors stay with in certain bounds. This approach will be taken 
here. The Jacobi fields are here treated outside the symmetric hyperbolic system following work 
in [23]. 

Since the reference solution is conformally flat its Weyl tensor and Cotton- York tensor vanish. 
However the second fundamental form and the Schouten tensor do not vanish and hence it is 
more difficult to obtain suitable bounds for the estimations that are to follow. Hence we use the 
second order Jacobi equation and replace the Schouten tensor with the Cotton- York tensor in it. 
The conformal Jacobi equation (|7b[) is linear. Thus it is sufficient to consider initial data with 



z(0) = 1, Tf(0) = e£, z(0) = Xkk- 
Furthermore, let X = z(0). We split the Jacobi fields into 

rf = TfoV^ + zm? 

where z 2 = —h(r),rj) and h{m,m) = — 1. Denoting d/dr by D, we can derive 

z = 1 > = {h ik + L ik )m m z, 



z 



where we have used the Cauchy-Schwartz inequality and the conformal Jacobi equation (|7bl 
We observe the identity D(L(rj, e k )) = Y(v, r), e k ) and its integral form 



o 



Yojkrfda = L(rj, e fc )(r) - L(r?, e fc )(0). 
This allows us to replace the Schouten tensor in the above inequality and obtain 

'z > E(m,m)(T)z(T) + [ (Y j k zm j )(a)m k (r)da + {zL ik m i )(0)m k (T). 



Observing that m is a space-like unit vector and employing a Cauchy-Schwartz type argument 
in combination with the existence and stability theorem ([1]), we can see that there exist bounds 
such that 

—K < £?(m,m)(r) < K — L < % jk m? (a)m k (T) < L. 
For the Schouten tensor on the initial surface we use 

-Q<L ik zm i (0)m k (T) < Q' , 

with Q, Q' > 0. This simplifies the inequality to give 



z > -Kz(t) - L ( z(cr)da - Q. 
Jo 
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One now follows the argument of |23j . by introducing a family of functions y e (r) satisfying 

y £ = -{K + e 2 )y £ (T) -L f y e (a)da - Q (30) 

Jo 

with initial data y e (0) = 1, y £ (0) = X. The idea is to show that z{t) cannot vanish before 2/ £ (t), 
in particular when e = 0. By estimating a lower bound for the time r y when y e= o vanishes, we 
get a lower bound for the time before the solution can develop conjugate points. 

We briefly recall the steps in the argument, whose details can be found in [231 IS]- We drop 
the subscript e. One defines R = z/y and W = zy — yz — y 2 R. From this one derives 

W{t) > L [ y(r)y(a)[R(T) - R(a)]da + Q V {t)[R{t) - 1]. 
Jo 

There must be an interval J = [0, T e ] such that W > and hence that R > 1 and W > hold 
on J. This implies that as long as y does not vanish, we can extend J further while R > 1 holds. 
Hence z > y > on J. By continuity in e, lim e ^o T £ — Tq > 0. Thus we only need to focus on 
y = y e =o, for which integrate ([30]) twice to get 

y{r) = \-XT-\Qr 1 -K f (r - a)y(a)da - L [ \{t - afy{a)da. 
1 Jo Jo 1 

Observing that y[r) < y(0) + Xt and using it to obtain upper bounds for the integrals we get 

1 1 LX 

y(r) > Y(r) =1-Xt--(Q + K)t 2 - -(XK + L)r 3 - — r 4 

For the reference solution all constants vanish so that z > y > 1 for all r, which agrees with 
z = 1 + t 2 /4. If one of the constants X, Q, K, L is non-zero, we get exactly one positive root, as 
all points of inflection have to lie in the second and third quadrant. 

Without loss of generality we can set all of them to R by appealing to Kato's theorem. Suppose 
we now fix t — T suitably beyond the values at future infinity of the space-time we would like to 
perturb. Then F(T) = 1 - A{T)R - B(T)R 2 , for some constants A(T), B(T). It is clear that we 
can always choose the perturbation of the initial data suitably small in Kato's theorem, so that 
the bound R guarantees that f(T) > and hence that the perturbed solution will not develop 
conjugate points before the chosen time r = T. 
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